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Abstract 

The semigroup decomposition formalism makes use of the functional model for C.q 
class contractive semigroups for the description of the time evolution of resonances. For a 
given scattering problem the formalism allows for the association of a definite Hilbert space 
state with a scattering resonance. This state defines a decomposition of matrix elements 
of the evolution into a term evolving according to a semigroup law and a background 
term. We discuss the case of multiple resonances and give a bound on the size of the 
background term. As an example we treat a simple problem of scattering from a square 
barrier potential on the half-line. 

1 Introduction 

Originally formulated for the analysis of scattering problems involving solution of hyperbolic 
wave equations in the exterior domain of compactly supported obstacles, the Lax-Phillips 
scattering theory 1 was developed as a tool most suitable for dealing with resonances in the 
scattering of electromagnetic or acoustic waves. Subsequent to its introduction by Lax and 
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Phillips, various authors have contributed to further development of the theory 2 ' 3,4 ' 5,6 . No- 
table recent additions were made by Sjostrand and Sworski 7 who extended the scope of the 
theory to include general classes of semibounded, compactly supported perturbations of the 
Laplacian in the wave equation, and by Kuzhell, via the development of a formalism pro- 
viding conditions for the application of the Lax-Phillips structure to an abstract form of the 
wave equation 8 and to certain classes of Schrodinger operators 9 . In addition, Kuzhell and 
Moskalyova 10 applied the Lax-Phillips theory in the analysis of scattering systems involving 
singular perturbations of the Laplacian. 

Several recent papers have dealt with the adaptation of the Lax-Phillips theory to quantum 
mechanical scattering problems. An early work in this direction is Ref. 11,12,13. A general 
formalism was developed in Ref. 14 and subsequently applied to several physical models in 
Ref. 15,16,17. Such efforts to adapt the Lax-Phillips formalism to the framework of quantum 
mechanics are motivated by certain appealing features of the Lax-Phillips theory. One of 
these features is the fact that the time evolution of resonances in this theory is given in terms 
of a continuous, one parameter, strongly contractive semigroup {Z(t)} t >o 

Z(ti)Z(t 2 ) = Z(*i + t 2 ), *i,t 2 >0. 

If TL is a (separable) Hilbert space corresponding to a particular scattering system and 
{U(t)} te iR is a unitary group defined on H describing the evolution of the system, the basic 
premises of the Lax-Phillips theory include the assumption of the existence of an incoming 
subspace T>^ and an outgoing subspace T> + with respect to {U(t)}i G iR which are assumed 
furthermore to be orthogonal to each other. Denoting by P and P+ respectively the projec- 
tions on the orthogonal complements of 2?_ and V + in H, and letting K, = H © (P_ © 
the Lax-Phillips semigroup {Z(t)}t>o defined by 

Z(t) = P+U(t)P_ = P K \J(t)P K , t>0, (1) 

annihilates T>± and maps K, into itself. The subspace K, contains the scattering resonances 
and the Lax-Phillips semigroup {Z(t)} t >o describes their time evolution. In the Lax-Phillips 
framework resonances are associated with pure states in the Hilbert space H. 

A basic difficulty encountered in the work on application of the Lax-Phillips theory in 
quantum mechanics originates from the fact that in this theory the continuous spectrum 
of the generator of evolution is required to be unbounded from below as well as from above. 
Hence a formalism utilizing the original structure of the theory, such as in Ref. 14, is not suit- 
able for application to large classes of scattering problems in quantum mechanics (except for 
limited types of problems, such as the Stark effect Hamiltonian 17 , or problems in a relativisti- 
cally covariant framework 15 ' 16 , which can be analyzed by direct mapping to the Lax-Phillips 
structure. The case of a Schrodinger equation with compactly supported potential may also 
be analyzed within the Lax-Phillips framework through the use of the invariance principle of 
wave operators 18 ). The subject of the present paper is a theoretical framework, termed the 
semigroup decomposition of resonance evolution, developed with the goal of overcoming such 
difficulties. Proposed by one of the authors (Y.S.) of the present article 19 ' 20 , this formalism 
makes use of the Sz.-Nagy-Foias theory of contraction operators and contractive semigroups 
on Hilbert space 21 which, from the mathematical point of view, is the fundamental theory 
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underlying the Lax-Phillips construction through the notion of model operators for Co class 
semigroups (see Section [2]). 

The presentation of the semigroup decomposition formalism in Ref. 20 is based on the 
following assumptions: 

(i) We are considering a scattering system consisting of a "free" unperturbed Hamiltonian 
Ho and a perturbed Hamiltonian H, both defined on a Hilbert space TL. 

(ii) ess suppcr ac (Ho) = ess suppcr ac (H) = R + . For simplicity it is assumed further that the 
multiplicity of the a.c. spectrum is one. 

(iii) The M0ller wave operators fi^ = fi (Ho,H) exist and are complete. 

(iv) The S-matrix in the energy representation (the spectral representation for Ho), denoted 
by S(-) has an extension to a meromorphic function S(-) in an open, simply connected, 
region ScC such that £nM is an open interval in K. The operator valued function 6*(-) 
is holomorphic in £ n C + and has a simple pole (we generalize to the case of multiple 
poles in Section [3] below) at a point z = fx G E H C~ and no other singularity in E, the 
closure of E. 

It is shown in Ref. 20 that there exists a dense set A C 7i ac (H) and a well defined state 
ipn G Tt ac (H) such that for any g G A and any / G 7i ac the properties (i)-(iv) above induce, 
for positive times, a decomposition of matrix elements of the evolution U(i) in the form 



In a sense to be made precise in the next section the second term on the right hand side of 
Eq. ([2]) originates from an evolution semigroup of Lax-Phillips type and the eigenvalue of the 
generator of this semigroup is exactly /i, i.e., the point of singularity of the S'-matrix. The 
quantity R(g, f;t) on the right hand side of Eq. ([2]) is what we shall call a background term. 
We note that if in Eq. (|2|) we choose / to be orthogonal to tpn then the exponentially decaying 
semigroup term (second term on the r.h.s. of Eq. ^fy ) vanishes. We call ip^ an approximate 
resonance state and note that the characterization of ipu, as an approximate resonance state 
rather than as an exact resonance state stems from the fact that one can show (see Ref. 20) 
that there is no choice of g and / that makes the backgound term R(g, /; t) vanish. 

An explicit expression for the approximate resonance state ip^ is provided in Ref. 20. 
It is shown there that, if we denote by {\E~)} E( z^+ the set of outgoing solutions of the 
Lippmann-Schwinger equation (using Dirac's notation), then ipu is given by 



Following the introduction of approximate resonance states, the present paper discusses some 
generalizations. Thus, in Section [3] we assume that the region E n C~ contains multiple 
resonance poles of the S'-matrix <S(-), say at z = fj,i,...,fj, n and obtain the form of the 
expression for the approximate resonance states and semigroup decomposition of evolution 
matrix elements in this case. In particular, we apply the semigroup decompostion to the 



{g, U(t)/) Wac(H ) = R(g, /; t) + a(g, /z)(Vv> /)w oc (H)e' 



t > 0. 



(2) 




(3) 
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survival amplitude, a central notion in the characterization of the time evolution of resonances. 
Theorem [5] below then provides an a priori upper bound on the size of the background term 
in this case. 

As a final remark we note that a modification of the Lax-Phillips theory was recently used 
by H. Baumgartel for the description of scattering resonances in certain quantum mechanical 
problems 22 (see also Ref. 23). In particular, the assumption of orthogonality of T>±, essential 
in the context of the original Lax-Phillips formalism, is replaced in Ref. 22 by the require- 
ment that an incoming subspace T>- and an outgoing subspace V + exist and the respective 
projections commute. The modified assumptions on D±, accompanied by certain assump- 
tions on S'-matrix analyticity properties, result in a modified Lax-Phillips structure which 
is then applied to the Friedrichs model, leading to the construction of appropriate Gamow 
type vectors 24 associated with scattering resonances. The framework presented in Ref. 22 
has several points of intersection with the semigroup decomposition formalism discussed in 
the present paper. The nature of these relationships will be discussed elsewhere. 

The rest of the paper is organized as follows: In Section [2] we describe the formalism 
providing the semigroup decomposition of resonance evolution starting with a short discussion 
of the functional model for Co continuous contractive semigroups followed by a description 
of the semigroup decomposition formalism introduced in Ref. 19,20. In Section [3] we extend 
the framework of Ref. 19,20 to the case of multiple resonances and, furthermore, find an 
estimate on the size of the background term in the expression for the time evolution of the 
survival probability of a resonance. In Section 0] we analyze a simple but illuminating example 
involving a one dimensional model of scattering from a square barrier potential. Section [5] 
contains a short summary of the contents of the paper and some indication on further possible 
courses of investigation. 



2 The semigroup decomposition for 
resonance evolution 

2.1 Classification of contractive semigroups 

Several distinct classes of contractive semigroups are identified within the framework of the 
Sz.-Nagy-Foias theory. Let {T(t)} t >o be a strongly contractive semigroup defined on a Hilbert 
space H. The classes Co., Co, C, C.\ are defined by 

{T(t)} tm+ E C Q . if T(t)h - 0, Vh E H 

{T(t)} tm+ E Co if T*(t)h -» 0, V/* G H 

{T(t)} tm+ E Ci. if T(t)h t^O, V/iG«, /i/0 

{T(t)} teM+ E Ci if T*(t)h A0,VkH,M0 

The classes C a p with a, (3 = 0, 1 are then defined by 

C a(3 = C a .nC.p, a, (3 = 0,1. 

The semigroup {Z(i)} tgK + describing the time evolution of resonances in the Lax-Phillips 
theory is readily characterized by the fact that {Z*(t)} fgK + belongs to the class Co • The 
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structure of the Lax-Phillips outgoing spectral (and translation) representation is then deter- 
mined by that of the functional model 21,25 for Co class semigroups provided by the Sz.-Nagy- 
Foias theory. We say an operator A is a model operator 25 for a given class C of operators 
if every operator in C is similar to a multiple of a part of A (a part of an operator A is a 
restriction of A to one of its invariant subspaces). By a functional model we mean that the 
model operator for a given class C has a canonical representation on suitable function spaces. 
For a Co class semigroup {T(t)}t>o the associated functional model is essentially obtained 
through a procedure of isometric dilation of the cogenerator of {T(t)} t >o and the similarity 
mapping to the functional model is in fact a unitary transformation. 

2.2 The functional model for C. semigroups 

We turn now to a brief description of the functional model for semigroups in the class Co 
. Denote by C + the upper half of the complex plane and let Hjj-(C + ) be the Hardy space 
of vector valued functions analytic in the upper half-plane and taking values in a separable 
Hilbert space N . The set of boundary values on R of functions in Hjj-(C + ), denoted below 
by Hjjj- + (R), is a Hilbert space isomorphic to H^-(C + ). In a similar manner the Hardy space 
of J\f valued functions analytic in the lower half-plane is denoted by Hjj-(C~) and Hjj-_{R) 
is the isomorphic Hilbert space consisting of boundary values on R of functions in Hjj-(C~). 
Define {u(t)} t eR, a family of unitary, multiplicative operators u(t) : L^-(R) i-> -Ljy-(R) by 

[u(t)f](a) = e-^fia), f G L 2 N (R), a G R. (4) 

Assume that {T(t)} t >o is a Co class semigroup defined on a Hilbert space /C. Let the 
semigroup {T(t)}t>o, defined on a Hilbert space fC, be the functional model for {T(t)}t>o and 
let W : K. i— > K. be the similarity transforming {T(t)} t >o into its functional model {T(*)}t>o 
i.e., T(t) = WT(i)W -1 . Then there exists a Hilbert space N such that JC is a closed subspace 
of -ffjy- + (R), W is unitary, and the functional model is given by 

f(t) = WT{t)W* = P t u*{t)\K, t>0. (5) 

Here is the orthogonal projection from iijy- + (R) onto /C, the subspace K, is given by 

K = 4 + (»)ee T (.)4 + (R), (6) 

and 6t(-) : Hjj- + (R) h-> H^-JR) is an inner function 26,29 ' 30 for H^ + (R) (depending, of 
course, on {T(t)} t >o) i.e., an operator valued function with the properties: 

1. For each a G R the operator ©x , (cr) : M ^ M is the boundary value at a of an operator 
valued function @t(-) analytic in the upper half-plane. 

2. \\O t (z)\\m < 1 for Imz > 0. 

3. 0t(c), a G R is, pointwise, a unitary operator on J\f. 
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The operator valued function Oy(-) is, in fact, the characteristic function 21 of the cogenerator 
of the semigroup {f(t)} t > (or {T(t)} t > ). 

Let .P+ be the orthogonal projection of L^-(R) on .£fj^ + (IR). The Toeplitz operator with 
symbol u(t) (see, for example, Ref. 26,27 and references therein), is an operator T n ( t ) : 
fljy + (K) h+ fl^ + (M) defined by 

T u(t) f= f P + u(t)f, feHfr + (R). (7) 

We note that {T u (^}t>o is a strongly contractive semigroup on fZjy-, (R) (see, for example, 
Ref. 1,19,21,28. Taking the conjugate of T(t) in H^ + (R) and using Eq. ([5]) one finds that 

f*(t) = WT*(t)W* = T u{t) \K, t>0. (8) 

It follows from the discussion above that the Lax-Phillips semigroup {Z(t)}t>o has a 
functional model in the form of Eq. ([H]) (recall that {Z*(t)} t >o is a Co class semigroup), i.e., 
if we denote the functional model for {Z(t)}t>o by {Z(t)}t>o then we have 

Z(t) = WZ(t)W* = T v(t) \K, t>0 (9) 
where K, C Hj^ + (M) is an invariant subspace for {T u (t)}t>o given by 

t = Hfr + (R) Q e z (.)iZfr+(R) (10) 

and the inner function Z (-) and the Hilbert space M are determined by {Z(t)} t >o- A 
semigroup {Z(t)} 4 >o of the form given by Eq. © and Eq. (fTU|) is referred to in Ref. 20 as a 
Lax-Phillips type semigroup. 

A central theorem of the Lax-Phillips theory, corresponding to an important result in the 
Sz.-Nagy-Foias theory relating the spectrum of a completely non-unitary (cnu) contraction to 
points of singularity of the characteristic function states the following 

Theorem 1 Denote by B the generator of a Lax-Phillips type semigroup {Z(t)}t>o- Iflm/i < 
0, then [i belongs to the point spectrum of B if and only ifQ* z (~fi) has a nontrivial null space. 

We note that the analytic continuation of @z(z) to the lower half-plane is given by 

Qz(z) = f (GHz))- 1 , lmz<0 

and so a null space for Q* z (Jl) implies the existence of a pole for &z( z ) at z = fx. In the case 
of the Lax-Phillips theory the characteristic function Qz(') for the Lax-Phillips semigroup is 
identical to the Lax-Phillips 5-matrix and its poles are the scattering resonances. As will be 
seen below, the situation is a bit more involved in the semigroup decomposition formalism. 

We do not elaborate here further on the relations between the functional model for C.q 
semigroups discussed above and the full structure of the Lax-Phillips spectral representations 
and wave operators. The reader is referred to Ref. 1,21. 
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2.3 The semigroup decomposition 

In order to apply the functional model for C.q semigroups , which is at the heart of the 
Lax-Phillips structure, to the description of resonance evolution it is necessary to relate, for 
t > 0, the evolution U(t) defined on the Hilbert space H of the scattering problem to the 
Toeplitz evolution semigroup T u (t) of Eq. (|7|) defined on H^- + (W) and then restrict the latter, 
according to Eq. ([9j), to a subspace K, of tij^, (R) associated with an appropriate inner 
function Q^(-). In the framework of the Lax-Phillips theory this relation is guaranteed by 
the special properties of the Lax-Phillips incoming and outgoing subspaces T>± (with 2?_ and 
V + denoting, respectively, the incoming and outgoing subspace), since in this case the Lax- 
Phillips semigroup is a C.q semigroup. However, for many quantum mechanical scattering 
problems one usually cannot find subspaces with the properties of T>±. A way of overcoming 
this difficulty, proposed in Ref. 19 is to combine the standard functional model for C.q 
semigroups with the notion of a quasi-affine mapping (see, for example, Ref. 21, Pg. 70): 

Definition 1 (Quasi-affine mapping) A quasi-affine map from a Hilbert space 7i\ into a 
Hilbert space TCq is a linear, one to one continuous mapping ofTCi into a dense linear manifold 
in TCq. If A £ B(TC\) and B £ B(TCq) then A is a quasi-affine transform of B if there is a 
quasi-affine map 6 : 7i\ i— > TCq such that OA = B#. 

The following theorem is proved in Ref. 19 for a scattering system consisting of unperturbed 
and perturbed Hamiltonians, respectively Ho and H, having semibounded continuous spec- 
trum: 

Theorem 2 (Outgoing/Incoming contractive nesting) Let Ho and H be self-adjoint 
operators on a Hilbert space 7i. Let {U(i)}t g R be the unitary evolution group on 7i generated 
by H [i.c, U(i) = exp(— iilt)]. Denote by TC ac (ilo) and Tl ac (H) , respectively, the absolutely 
continuous subspaces of Hq and H. Assume that the absolutely continuous spectrum of Ho 
and H has multiplicity one and that ess Supp er ac (Ho) = ess Supp o~ ac (H) = M + . Assume 
furthermore that the M0ller wave operators fi ± = fi^Ho,^ : W ac (Ho) i— » 7i ac (H.) exist and 
are complete. Then there are mappings Cl± : H ac (H) i — ^ Hl(R) such that 

(i) Cl± are contractive quasi-affine mappings ofH. ac (H) into H+(M). 

(ii) For every t > the evolution U(t) is a quasi-affine transform of the Toeplitz operator 

T n ( t ) via the mapping Cl± i.e., for every f £ W ac (H) we have 

fl±U(t)f = T <t) fl±f t>0. (11) 

□ 

We call the triplet (7i oc (H), H^_(M), 0_) the incoming contractive nesting of W ac (H) into 
H^_(M.) and denote fi n = &-f. Similarly, the triplet (W ac (H), H+(M), Q + ) is the outgoing 
contractive nesting of H ac (H) into iJ? (R) and we denote f ou t = &+f- 
Define 

3, d MCl* + Hl(R). 
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Then, since fi+ is quasi-affine, the linear space C TC ac (H) is dense in 7i oc (H). Moreover, 

since Cl+ is one to one, for each g £ there is a unique 5 £ H^.(R) such that g = &* + g. 
We note that in Ref. 20 a dense set , analogous to , is defined somewhat differently, 

i.e., = Cl^_Cl + Ti ac (ii) . However, it will be seen below that the definition of above, 
unlike that of A^ + , allows for a full characterization of approximate resonance states. Using 
Theorem [2] we have, for every g £ and / £ H ac (H) and for t > 

(g,U(t)f) Hac(u) = (n;&u(t)/) Woc(H ) = 

= (9,T u (t)£l + f) H 2Qgr) = (g,T u (t)fout)HlQSL)i ( 12 ) 

Following the definitions of the incoming and outgoing nestings of H ac (H) into JET£(R) it 
is natural to define the nested S-matrix 

Snest = ^+^_ • 

Let U : TC ac (Ho) 1— > L 2 (R + ) be the unitary transformation of TC ac (Ho) onto the spectral 
representation for Ho (also called the energy representation for Ho). If S = (fi - )*fi + is the 
scattering operator associated with Ho and H then §(■) : L 2 (R + ) 1— ► L 2 (R + ) defined by 

S(-) = USU* (13) 

is the energy representation of the S-matrix. Let : L 2 (R) 1— ► L 2 (R) be the orthogonal 
projection in L 2 (R) on the subspace of functions supported on R + and define the inclusion 
map I : L 2 (R+) \-> L 2 (R) by 



(//)(-)= \r h (14) 




Then the inverse I 1 : Pjj+L 2 (M) 1— > L 2 (R + ) is, of course, one to one on P R +L 2 (R). Let 
9 : Hl(R) 1 — ► L 2 (R+) be a map given by 

9f = L~ 1 P R+ f, f€Hl(R). (15) 

By a theorem of Van Winter , 9 is a quasi-affine transform mapping H\{R) into L 2 (R+). 
The adjoint map 9* : L 2 (R + ) ^ H+(M) is then also a contractive quasi-affine map. An 
explicit expression for 9* is provided by the following lemma 19 : 

Lemma 1 Let I : L 2 (R + ) L 2 (R) be the inclusion map defined in Eq. |i^[ ). Let P + be the 
orthogonal projection o/L 2 (R) onto ii+(R). Then for every f £ L 2 (R+) we have 

9*f = P+If, f £ L 2 (R+) . (16) 

□ 
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It is shown in Ref. 19 that the nested S-matrix can be expressed in the form 

S nest = 9*S(-)(8*)- 1 . (17) 

Following Ref. 20 we now use assumption (iv) in Section [U The S-matrix S(-) is then the 
restriction of its extension <S(-) on R + . Under these assumptions <S(-) has, in the region E, a 
representation of the form (see Ref. 20) 

S(z) = B^(z)S'(z), z e £ 

where _ 

B„(*) = ^, zeC\{ M } (18) 

Z — jU 

and <S'(-) is analytic and has no zeros in S. Restricting <S(-) to the positive real axis we obtain 

S(E) = B IX (E)S'(E), E>0 (19) 

where by definition B^E) d = B^(E) and S'(E) d = S'(E) for E > 0. We note that both S'(-) 
and B„(') are considered here as multiplicative unitary operators on L 2 (R + ) (moreover, they 
are pointwise unitary a.e. for E > 0). Moreover, B^(-) can be regarded as a multiplicative 
operator on L 2 (R). In fact, considered as a multiplicative operator on H 2 _(R) C L 2 (R), 
is a Blaschke factor (the definition of Blaschke products and Blaschke factors can be found, 
for example in Ref. 29,30, see e.g., Eq. (I29p below). Such a factor is the simplest example 
of an inner function for H 2 _(R). We make use of this fact through the following proposition, 
not stated as such, but implicitly used in Ref. 20: 

Proposition 1 Let B^(-) : L 2 (R + ) i-> L 2 (R+) be defined by B^E) = B^{E), E > where 
BJ-) is defined in Eq. Wto. Let 9* : L 2 (R+) i-> Hl(R) be the adjoint of the map 9 defined 



in Eq. [75j) . Let C H+(M) and Kp C H 2 ^ (R) 6e subspaces defined by 

= f H 2 + (R) Q B^HKR), tp = f H 2 (R) Q Bp(-)H 2 _(R) , 

where B-p(z) = (z — fj){z — 7i) _1 and denote by Pg and Pg-_ the orthogonal projections of 
L 2 (R) on K,^ and Kp respectively. For every f G L 2 (R + ) we then have 

9*B^f = B^9*f + P t B»P t _tf . (20) 

here 9* : L 2 (R + ) i-> F 2 (R) and 9* f = P^Lf with L defined in Eq. and P_ the orthogonal 
projection of L 2 (R) onto H 2 (R). □ 

Proof: Using Eq. (|16p in Lemma [T] we get 

0*B„f = P+LB^f = P+B^Lf = P + B^(P + + Pj)If = P + B^9*f + P+B^Tf 

Eq. (|20p then follows from the fact, proved in Ref. 20, that P+B^P~ = Pg- B^P^_ and from 
the property of B^(-) of being an inner function for iT? (R). 
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We note that since B u (-) is an inner function Eq. ([9]), (I10p and Theorem [T] imply that 

T <t)P tli = &®P U = e**P^ t>0. (21) 
Combining Eq. (JTTJ) , Eq. (PJ) and Eq. (J20|) we obtain 

/out — driest Jin — V ) Jin — V -£>>-> \y ) Jin — 

= Ba9*S'{9*)~ 1 f in + P t B.P^Ts'ie^fin • (22) 



Using the decomposition of f ou t from Eq. ([22l in the r.h.s. of Eq. (fT2|) and applying Eq. (|2Tj) 
we obtain the semigroup decomposition for t > of the time evolution corrsponding to the 
resonance at z = fJ> 



(flS U(£)/) Woc ( H ) = (g,T u ( t )fout) H l(R) = 

= (g,n(t)^r5'(r)-V m )m (R) + e -^(g,^r5'(r)-V m )m m - (23) 



As is seen above in Eq. (|2ip . the exponential decay in the second term on the r.h.s of Eq. (|23[) 
originates with the semigroup Z(i). The first term on the r.h.s. of Eq. (|23p is the background 
term and is responsible for deviations from a purely exponential decay law. 

3 Approximate resonance states 

It is an interesting fact that the semigroup decomposition described in the previous section 
associates a unique state in TC ac (H.) with a resonance pole at z = n (Im fj, < 0). The following 
theorem is proved in Ref. 20 

Theorem 3 (approximate resonance state) Under the assumptions of Theorem [IJ let 
S : L 2 (R + ) i — ^ L 2 (1R + ) be the S-matrix in the energy representation defined in Eq. $13\) . 
Assume that S(-) is the restriction to R + of a function S(-) meromorphic in an open region 
£ with a single, simple pole at a point z = fi, /i G S D C - . For any f £ TL ac (H.) define 
fout = £l+f and f n = 0_/. There exists a unique state ip n E H. ac (H) such that 

f out = S^rS^T 1 /™ + /)w ac (H) ^ (24) 

where 0* is the map given by lemmaUl B u is given in Eq. U8\) . S' is defined by Eq. U9\) and 
x^ £ fl+(R) is given by x u (a) = (a — /u) -1 , ueK. □ 

Combining Eq. (I24p and Eq. ([23]) we can write the semigroup decomposition in the form 



( 5 ,U(t)/) Wac(H ) = (g,u(t)B lt 0*S\9*)- 1 fin) H * 

+ l -^e-^ t (g,x fl ) Hlm (^,f)n ac (H), t>0 (25) 



7T 
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where g G and g = (O^.) The eigenstate rc^ G iT 2 (R) of the semigroup Z(t) providing 
the exponential decay of the second term on the r.h.s. of Eq. (I25D is called below the Hardy 
space resonance state. The state W e W ac (H) whose existence is implied by Theorem [3] is 
called approximate resonance state. We observe that if in Eq. (|25|) we choose / G W ac (H) 
orthogonal to then the second term on the r.h.s. of that equation is identically zero. 

Denote by {|-E~)}_egr+ the set of outgoing solutions of the Lippmann-Schwinger equation. 
For every / G 7i oc (H) we have 

(U(n-)*f)(E) = (E-\f) 1 E G R + (26) 

It is shown in Ref. 20 that an explicit expression for the approximate resonance state ifjn is 
given by 

% = I dE-^—\E-). (27) 
2m J R + E - ii 

In this section we explore several properties of approximate resonance states ip^. Our first 
step is to extend the discussion above to the case of multiple resonances: 

Theorem 4 (multiple resonance case) Under the assumptions of Theorem^ let S : L 2 (M" t 
L 2 (R + ) be the S-matrix in the energy representation defined in Eq. I113\) . Assume that S(-) is 
the restriction to R + of a function <S(-) meromorphic in the open region S with n simple poles 
at points z = m, i = 1, . . . , n, \ii G E Pi C _ . Then there exist n distinct states {V^}i=i,...,n, 
V>u G H oc (H), smc/i £/tcrf /or every / G W ac (H) iue /iaue 



fout — ^ti...u„ 



3=1 

where 



*S>(9*r'f in + E ^ II . /)«-CB) (28) 



n 

B ^) = nF?' (29) 

In .Eg. j&S)) 5"(-) is i/ie restriction to M + of a function S'(-) analytic in S and having no poles 
in S and x^. (cr) = (<r — fJ-j)^ 1 , o~ G R. TTie states ^ , j = 1, . . . ,n are given by 

€ 3 = I ^n r ft / i^">- (3o) 

□ 

Proof: Assume that <S(-), the extension of S'(-) from M + into SUR + has n simple poles in 
£nC~. Then, applying the same arguments as in Ref. 20, we find that <S(-) can be factorized 
in S in the form 

S(z) = £ Ml ... M „(z)S'(z) 
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where ^ Atl ...^ n , defined in Eq. (f29l) . is a finite Blaschke product and S'(-) has no poles in S. 
In addition we have, of course 

S(E)=B^ n (E)S'(E), E>0. 

The semigroup decomposition then follows exactly as in Section [T] with S Ml ... /Jn and -B^...^ 
replacing B^ and respectively. For the resonance term in Eq. (|22[) we get in this case 

Recalling that 

1 f°° 1 

(Pf/) (*) = 5-7 / A/ 7— rr/(o-'), / G L 

z7Ti a — a' + 



2;lro ), CJG 



and 

_ 1 poo 1 

<"><*> = Si/. iE E^Tm ! ^ /£i2 < R+) 
(see Ref. 19) we obtain 

{P+B tXl ... lin PJFS'{0*)- 1 f in ){o) = 

— 1 f 00 f 00 1 J-L-rr'—Ti 1 

= t4/ / d*' 1 TTV-^^ -, S'{E){(6*y l f m ){E) 

^\iu^ 2 [^_ r dE f[ V-^ 1 S>(EW*T 1 U(E) = 

±3 it a- Hj J a± ^ -mE- m 



3=1 



T7 /"OO 



It It i I >■ nQQ -i 'fc t— i 

j=l i=i ^ •'O O i=i 

71 I T i /I n OO -« 71 T— • 



i=i 



I Ij It i /'OO 1 7— ■ 

e^mr^/ ^nfff< E -i/> 



where x flj G -H+(M) is the Hardy space resonance state corresponding to //j i.e., x^Aa) = 
(a — ^j)^ 1 - Defining the states ipjj., j = 1, . . . n according to Eq. ([30]) we obtain 

n ,j I n 

(p + w ..^p-r^(r)- 1 /«)(a) = E^ Ilir^r^'^w^W ( 31 ) 
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This proves Theorem 01 



We observe that Eq. (|30p is a generalization of Eq. (|27p . Hence tp^. is the approximate 
resonance state corresponding to the pole of <S(-) at z = fij. Combining Eq. (j3Tj) and Eq. 
(|23p we get the semigroup decomposition for the multi-resonance case 



(5,U(i)/) Woc(H) = (9,u(t)B (i ,.. )1 „r5'(r)- 1 / m ) ffl(E ) 

j=l i=l J 

The approximate resonance states in Eq. (I30p and semigroup decomposition of Eq. (|28[) 
and Eq. ([32]) depend, of course, on the region S. If {Mj}j=i,...,n are the poles in SnC" of the 
meromorphic extension S(-) of the 5-matrix S(-), the approximate resonance state defined 
in Eq. (|30p f° r a resonance at z = fij is therefore denoted by tp^.. However, for certain 
arguments the exact form of £ is irrelevant and it is useful to define the notion of an n'th 
order approximate resonance state: 

Definition 2 (n'th order approximate resonance state) If the number of poles ofS(-) 
entering into the definition of the approximate resonance state tp^. in Eq. \30\) , not including 
fij itself is n we say that ipj}. is an n 'th order approximate resonance state for the resonance 
at z = fij. In particular, regardless of the exact nature of the region S, the zero'th order 



approximate resonance state is always defined to be given by Eq. \21) with fi = fij and is 
denoted by . 

Remark: Note that in general there are many choices of the n resonance poles (different 
than fij) included in the construction of what we call an n'th order approximation ip^/ . In 
cases that the nature of the region S is irrelevant and only the order of the approximate 
resonance state is significant we replace the notation tp^ by ipjf^ , where n is the order of the 
approximate resonance state considered. 

The semigroup decomposition and approximate resonance states for the multi-resonance 
case possess some interesting properties. For example, we have 

(<,^ fc W(H) = (< ) ^£ ) W(H)= / dE- ' ' 
and, in particular 



E — fij E — fi k 



W^l IIL.CH) = U«~(H3 = ll«~CH> = L dE ■ (33) 

JM.+ 1-° H'j\ 

We see that, although the definition of ift%. in Eq. ([3"U]) depends on all of the poles {fij}j=i t ... t n C 
£ n C _ , the scalar product of ipj^. and ip^ depends only on fii and fij. In fact, if <S(-) can 
be extended to a meromorphic function in a region £' D £ (we keep the notation S(-) for 
the extended function) and <S(-) has now m > n simple poles in £' n C _ we may calculate 
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approximate resonance states of order m — 1 for all resonances in £' according to Eq. (I30p . 
However, for G £ we would still have (^', ^)K ac (H) = (^>^ fc )w oc (H) i- e -> scalar 

products (and norms) are independent of the order of the approximate states when we enlarge 
the region E. In particular we have ||^ ac (H) = IIV^ \\H ac {H) f° r every region £ containing 
Mi- 

An interesting question is whether the peculiar properties of scalar products and norms 
of the approximate resonance states mentioned above characterize also the time evolution of 
these states. We shall see below that, at least partially, the answer to this question is positive. 
For this we consider one of the basic notions associated with resonance evolution, i.e., that of 
the survival amplitude 

z , \ def Wac(H) /nA\ 

Making use of Eq. (|3U[) and (|33|) we get a simple expression for this quantity 
A (o) (t) = W^WZ 2 m) [ lr , 1 l9 e~ im , t > 

where || is given in Eq. (1331) . Again, we see that the expression for the survival amplitude 
for the approximate resonance state tp^. depends only on the pole at z = fij and has the same 
form as for a single resonance. This suggests that the semigroup decompostion of the survival 
amplitude for the multiple resonance case is similar to that of a single resonance. When 
combined with an important characterization of approximate resonance states in the form of 
Lemma [2] below, such considerations lead to the following useful a priori estimate on the size 
of the background term in the semigroup decomposition of the survival amplitude: 

Theorem 5 Let A^y. (t), j = 1, . . . ,n be the survival amplitude defined in Eq. (34\ ) and let 
the background term . (t) be defined by the relation 



A^ (t) = R„ (t) + e-*«* , t > . (35) 



Then we have 



I II 4 



1/2 



\R,M< \ „ /( o)„4 - 1 ! • ( 36 ) 

where x^Aa) = {cr — jjLj)^ 1 is the Hardy space resonance state andtp^ E W ac (H) is the zero'th 
order approximate resonance state corresponding to the resonance at z = fij. □ 

Proof: We first have 

Proposition 2 For j = 1, . . . , n, let iftj}. be defined by Eq. A3(J\) and let A^ (t) be the survival 
amplitude defined in Eq. ffify . Then 

A^.(t) = U^^x^u^B^ t > 0. (37) 



14 



where x llj £ fl^.(K) is i/ie Hardy space resonance state and ip^) E H c(H) is the zero'th order 
approximate resonance state corresponding to the pole at fij, and S^.(-) is defined as in Eq. 
i.e., 

~ S 'n W = fz§7 = % (E 1 ) 5(E) . 

□ 

Proof of Proposition [2} We need first the following easily proved, but important, 
lemma 



Lemma 2 For j = 1, . . . , n, let ipjj. be defined by Eq. [30]) . Define 

n 

&ni...fi k ...n„(z) = Y\ 
1=1 



i=l ^ 



where {fJ-j}j=i,...,n are the poles of S(-) in £, and let x^. (a) = (<r — 1 . TTien we /iave 

In particular ip^ 1 ) = fl+x^ . □ 

Proof of Lemma [2} It is proved in Ref. 19 that, if fl^ are the M0ller wave operators, 
U : 7i ac (ilo) l— > L 2 (M + ) the mapping to the energy representation for Ho (see Eq. (|13p 
above) and 0* the map given in Lemma [l] then the quasi-afhne nesting maps fi± are given by 
Q± = 0*[/(f2 =F )*, hence we have f2+ = Q~U*6. Furthermore, by the definition of we have 

n F - 77 1 

(JW'ttl^^FTFT- SGM+ - (38) 

Moreover, according to Eq. (|26j) for every 6 L 2 (R + ) we have 

n-C7* ff = f dE\E-)g(E). (39) 

Applying Eq. ([39]) with g = &Bn 1 ...fe„. lin x lij and comparing with Eq. ([30]) proves the lemma 

■ 

Note that by Lemma[2]we have {p.*^)~ l ip^. = B )J , 1 ^p,._ lXn x li .. Hence by Eq. ([32]) we get 
A^,(*) = II<II«L ( h)(< 5 U(*)<)^c(h) = 



A |Im/i fc | " ^-7J. (^,^.)H ae (H) < t 

+ - 11 ... _ ... IL/.K 112 l^i-Mi-Mn^.^i^tlR) 6 

fc=l i=l^ fe ^ H^HWac(H) 



(40) 
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In the second term on the r.h.s. of Eq. (|40p we first separate the term with k = j and get 



-in k t 



E |Im n k \ ^ ft-ft (V^ fc , ^% )-H ac (H) 
fc=l ff i=l^ fe ^ ll^llw ac (H) 

~ 1^ - 11 „._„. 112 l^...Ai---M»^>%Jl^(R) e + e 

fe=l ^ i=l //fc ^ H^illWac(H) 

, t > 

Here, use has been made of Eq. (|42p below. The above expression can be further simplified 
since x^ k S /C^ fc = H+(M) Q £>^ fc (-).£f+(R) implies that for k ^ j we have 

(Bfil...fi j ...H n x fij •! x fi k )l{2 (R) = (BftkBfM 1 ...fij...fi k ...ti n x fi. j , x fj. k )l{2(M.) = 

where 

n 

def TT Z — 



Bfi 1 ...p,j...p, k ...fi n (z) — Y\ 
i=i 



i=l 

The first term on the r.h.s. of Eq. (|40p can also be simplified. We have 

(^,.. Aj ,..^x Mj ,n(t)^ 1 ...^r5'(r)-VJ, m )^(R) = 

= («(-t)B Ml .../i i ..., 1 „a; Mj ,B w ..., ln e*S'(e*)~VS,in)fla(R; 

= (^w---Ai---Mn w (-*)^' S W---^^*^'(^*) - VJ > m)H2( M ; 



(x w ,u(*)B w (9*5'((9*)-VS,iJfl«(») > * > 



Moreover, using Lemma [2] we find that 



Recalling that Eq. (j33j) implies that IIV^- ll"H K c(H) = llVyj llw ac (H) the proof of Proposition [2] 
is complete. ■ 

From Proposition [2] we see that, independent of the region E, the semigroup decomposition 
of the survival amplitude depends only on the zero'th order approximate resonance state. 
Comparison of Eq. ([371) and Eq. (f35l) gives 
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This expression for (t) is identical to the zero'th order background term we would get 

from Eq. (gSJ with f = g = . We now exploit this fact to obtain the desired estimate in 

Theorem Applying Theorem to the zero'th order approximate resonance state tp^ we 
obtain 



Now, since both Cl + and are contractive we note that Lemma[2]implies that , out I Iff 2 1 
||fi+^+ai/^ Hij-a (K) — W-Efij \ \h 2 (R) - I n addition in H+(M) we have x^ lj _L B^iJ^R). Therefore, 



liK ft S*' : l- 1 ,/, (0) II 2 I l Im ^| 2 ||,/.(0)||4 II,. ||2 || ||'J 

It is easy to verify that 



7T 



Mkm=/ ^-p^ = ^' ( 42 ) 



hence the inequality above can be written in the form 

//}^-l,/,(0). ||2 < (» 1,4 _ _|L/,(0)||4 \ _ ||-2 



^/*^.(r)->^ n ||^ (R) < M|X,||^ (R) - ||<1^ c(H ) ) ■ ( 43 ) 



Applying the Schwartz inequality to the r.h.s. of Eq. (|41l) and using the bound from Eq. (jl3 
we get the estimate in Eq. ([36]) . I 



As mentioned above the background term cannot be identically zero. Hence deviations 
from exponential decay of the survival probability are to be expected. In fact, it is easy to 
verify that the survival probability behaves for short times as |/4^, M .(i)| 2 = 1 — 0(t 2 ). Note 
that Eq. (f55|) implies that at t = we must have R^ (0) = 0. This is also seen from Eq. (|4T]) . 
since for t = we have u(0) = 1 and x^. _L 0^. (-)iT^.(R). Deviations from exponential decay 
are then due to the fact that x flj JL u(t)Bp.H+(M) for t > 0. 

4 Example: Scattering from square barrier potential 

In this section we apply the results of the previous two sections to a simple one dimensional 
model with a square barrier potential. Although simple, this model provides a good illustra- 
tion for the various results obtained above. In particular, we present numerical calculations of 
approximate resonance states of various orders accompanied with plots of the time evolution 
of the corresponding survival amplitudes and estimates of the size of the background term 
following from Theorem [5j 

The model we consider is a Schrodinger equation in one spatial dimension on the half- 
line M + with a square barrier potential. Thus we consider the free Hamiltonian Ho = — d% 
acting on L 2 (R + ) (where Hq is defined as a self-adjoint extension to L 2 (M + ) from the original 
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domain of definition D(— <9 2 ) = {4>{x) £ W| I = 0}) an d the full Hamiltonian is 

given by H = Hq + V where V is a multiplicative operator (Vf)(x) = V(x)f(x) with 



V(x) 



0, < x < a 
Vo, a < x < b 
0, b<x, 



where b > a > and we take Vo > 0. In this case there are no bound state solutions of the 
eigenvalue problem for H and we have <r(H) = er ac (H) = M + . In order to find the scattering 
states, calculate the S'-matrix and finally the approximate resonance states for this problem 
one solves the eigenvalue problem 

-dl^ E {x) + V(x)ip E (x) = Eil> E (x), EgR+ 

for the continuous spectrum generalized eigenfunctions iPe(x)- Imposing boundary conditions 
one finds that 

a±(k) sin kx, 0<x<a 
iPe(x) = < a 2 (k)e ik ' x + /32(k)e~ ik ' x , a < x < b (44) 
a 3 (k)e ikx + p 3 (k)e~ ikx , b<x 

where k = £ 1/2 and k' = y/E- V for E > V > or k' = iy/% - E for V > E > 0. The 
coefficients in Eq. (j44p are given by 32 



a 2 (k) 

/%(*) 
a 3 (k) 



1 

— ( 

2 

1 

2 & ' 



-ik'a 



k 

sin ka -\ — — cos ka a±(k) 
ik' 



ik a 



sin ka — cos ka aAk) 

ik' 



4 6 



-ikb 



[l + k'/k)e 



ik'(b-a) 



(sin ka H — — cos ka) 
ik' 



(45) 



k 



+ (1 - k'/k)e- ik,( - b - a \smka cos ka) 

ik 



ai(k) 



1 

4 6 



ikb 



k 



(1 - k'/k)e ik ' (b - a) (sin ka + -- cos ka) 

ik' 



(l + k'/k)e 



-ik'(b-a) 



(sin ka 



k 

Ik' 



cos ka) 



a\{k) 



with ai(k) to be determined by normalization conditions (see below). 

Given the full set of solutions {^Pe(x)}e£R+ f° r t ne continuous spectrum it is easy to 
find the sets {V'|K ::c )}egr+ °f solutions of the Lippmann-Schwinger equation corresponding 
to incoming and outgoing asymptotic conditions. Using Dirac's notation we have 32 



(x\E + ) 
(x\E~) 



-1 IpEjx) 

2% f3 3 (k) 
2i a 3 (k) 



(46) 
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Figure 2: Probability densities of approximate 
resonance states ijyjy , J = 1,2,3 for a = 2, b = 3, 
V Q = 10. Dashed lines - IV^O)] 2 , j = 1,2,3. 
Solid lines - |^ 9) (a:)| 2 , j = 1,2,3. 



Figure 3: Probability densities of the approximate 
resonance states for fi' 3 = 17.4652 — i4.4029. Dot 

dashed line - \tf>^°) ( x )\ 2 > dashed line - (a;)| 2 , solid 
line - \ip { ' 3 } {x)\ 2 . 



The normalization conditions for the Lippmann-Schwinger states in Eq. (|46p determines 
a±(k) with the result a\(k) = (2nk)~ l l 2 . In the energy representation the S'-matrix is then 
given by 32 

x a 3 (k) _ 1/2 

s{E) -~Mk-y k ~ E 

We now have the ingredients for the calculation of the scattering resonances and the corre- 
sponding approximate resonance states. Note first that 03 (k) and (5^{k) can be extended to 
analytic functions in the complex k plane and as a result the poles of the analytic continuation 
of S(E) to the lower half-plane (i.e., across the square root cut along the positive real axis) 
are identified with zeros of the function (3^{k). For a resonance at a point z = fij in the lower 
half-plane below the positive real axis we set jij = E^. —iT^/2, with E^. being the resonance 
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energy and T llj its width. 

Using Eq. ()27[) and the expression for the outgoing Lippmann-Schwinger eigenfunctions 
(x\E~), Eq. (fill) . (j4"B"j) and (|46p . the zero'th order approximate resonance states for the square 
barrier problem can be calculated numerically. Considering a larger number of resonance poles 
we are able to calculate higher order approximate resonance states using Eq. (I30p . As an 
example we consider the three lowest energy resonance poles for barrier parameters a = 2, b = 
3, V Q = 10. These poles are located at \i x = 1.8213 - i0.0023, /i 2 = 7.0237 - i0.0564 and ^ 3 = 
14.2336 — iO. 8923. The zero'th order probability densities \ip^ (x)\ 2 , j = 1,2,3 for these poles 
are shown as dashed lines in Fig. [2] while the solid lines on the same figure correspond to the 
9'th order probability densities \ip^(x)\ 2 , j = 1,2,3, where the ten lowest energy resonances 
/jjj, j = 1, . . . , 10 are taken into account in Eq. (|30p . We observe the significant change in the 
probability density profile between the zero'th and 9'th order approximate resonance states 
for the resonance fi^, whose energy is higher then the barrier's energy, while the lower two 
states are essentially unchanged. Numerical calculations show that approximate resonance 
states converge in I? norm to a limiting state as a function of the order of approximation. 
An example is provided in Fig. [3] which shows the probability density \ip , (x)\ 2 for the 
resonance /ig = 17.4652 — M.4029, at barrier parameters a = 2, b = 2.1, Vq = 10, and for the 
orders n = 0, 8, 9. At present a rigorous criterion for the rate of convergence of approximate 
resonance states as a function of order is not yet established. 

Turning to a consideration of the time evolution of survival probabilities for resonances of 
the square barrier model, we first recall the fact that the time evolution of the survival prob- 
abilities of higher order approximate resonance states corresponding to the same resonance 
pole is independent of the order and is, in fact, identical to that of the zero'th order state. 
Bearing this in mind we may omit in our notation any indication of the region £ or the order 
n and set A^^(t) = A^s. (t) = A ,(n)(i). The time dependence of the survival probablility 

l^-i/Vi f° r ^ ne states coresponding to the lower resonance in Fig. [2] is shown as a solid line 
in Fig. [H The time evolution of | A/> Ml (^)| 2 follows closely an exact exponential decay law with 
a decay constant r w = 2|9f(^xi)|. This behaviour is reflected in the bound |i? Ml (t)| < 0.028 
on the size of the background term calculated using Theorem [5] The time development of 
\Ajjj (t)\ 2 deviates from the exponential law at a very short time scale, as is clearly seen in 
the insert in Fig. [H The behaviour of the survival probability for the other resonances in 
Fig. [2] (not shown in Fig. [T|) is similar. The short time deviations from exponential decay are 
related to the known Zeno effect. 

The nearly exact exponential decay law of the survival probability |yl^ Mi (t)| 2 is to be 
contrasted with the time development of \A^ , (t)\ 2 for the states in Fig. [3j The survival 

probability \A^ , (t)\ 2 is described by the dashed line in Fig. [TJ Deviations from an exponential 
decay law in this case are evidently larger. This conforms with the results of Theorem [5] which 
produces the larger bound \R^ (t)\ < 0.422. 
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5 Summary 



The semigroup decomposition formalism makes use of the fundamental mathematical theory 
underlying the structure of the Lax-Phillips scattering theory, i.e., the functional model for 
Co contractive semigroups, for the description of the time evolution of resonances. If the 
S'-matrix is meromorphic in a region £ and is known to have resonance poles there at points 
z = m, . . . , fi n S X, the semigroup formalism allows for the association of a unique Hilbert 
space state ip^.(x) G H ac (H.), j = l,...n with each resonance. The states ip^.(x), called 
approximate resonance states, define the decomposition of matrix elements of the evolution 
and are associated with its semigroup part. Theorem provides an upper bound on the 
size of the remaining background term. Depending on one's knowledge of the location of the 
resonance poles it is possible to calculate approximate resonance states of different orders. 
Numerical calculations show that the sequence of approximate resonance states appear to 
converge in 1? norm to a limiting function as a function of the order. However, rigorous 
criteria for the rate of convergence are needed. Another possible course of further investigation 
involves the study of relations between known frameworks for the treatment of the problem 
of resonances, such as the rigged Hilbert space method and the use of dilation analyticity and 
the formalism discussed in the present paper. 
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